Meraorandim M-2569 Page 1 of 1? 

Digital Computer Laboratory- 
Massachusetts Institute of Technology 
Cambridge 39 t Massachusetts 



SUBJECTS BIWEEKLY REPORT. DECEMBER 13. 1953 

Toll Jay ¥, Forrester 

From; Scientific and Engineering Computation Group 

lo MATHEMATICS. CODIUG AM) APPLICATIONS 

1»1 Introduction 

During the period covered by this report 238 coded programs were run on the 
time allocated to the Scientific and Engineering Computation (S&EC) Group. These 
programs represent part of the work that has been carried on in 28 of the prob- 
l^as that have been accepted by the S&EC Group, Progress on 23 of these problems 
is given below in teims of programming hours, minutes of computer time, and pro- 
gress reports as submitted by the programmers in question. 

Two new problems were initiated during this period. Problem #153 seeks to 
detemiine the gust response of a flexible swept-wing airplane under various 
conditions* This work is being carried out by K, Foss of the Aeronautical 
Engineering Department with the assistance of D, Sternlight of the Mathematics 
Department. In problem #157, Dr. A. Meckler of the Solid State Physics Group 
has developed i^outines for multiplying rectangular matrices of arbitrary order 
(up to 32)0 

A revised version of the programmed arithmetic (PA) section of the compre- 
hensive system (CS) has been developed to make available a significant amount 
of mistake anticipation, an unlimited number of cycle counters, and an extended 
and improved use of buffer registers. The new version of the PA will not be 
put into general use for at least one month but will be made available to 
programmers upon special request and under staff supervision. This time delay 
will provide for more extensive testing of the routines and for a more gradual 
transition from the old PA to the new* 

Dr. J.H.Laning Jr. of the Instnjmentation Laboratory spoke at the Seminar 
on Computing Machine Methods on December 1, 1953* He described a program (that 
he and I. Zierler have developed) capable of accepting mathematical equations 
expressed on Flexowriter punched paper tape in ordinary mathematical notation 
(within certain limits imposed by the Flexowriter) a-s input and automatically 
providing the desired solutions. 

The GS introductory programming course was begun on December 7» 1953 as 
announced. A report on the coarse is given below under problem #100« 
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1,2 Programs and Computer Operation 

The following summary is included as a guide for interpreting the abbreviations 
used below^ A more detailed description of the terms involved can be found in 
M-2497. 

a. The upper case letter following the problem number has the following 
significance: 

A implies the problem is NOT for academic credit, is UNsponsored* 

B implies the problem IS for academic credit, is UNsponsored» 

C implies the problem is NOT for academic credit, IS sponsored, 

D implies the problem IS for academic credit, IS sponsored. 

The absence of a letter indicates that it is an internal SSEC problem, 

b. DIG denotes the Division of Industrial Cooperation, 
DCL denotes the Digital Computer Laboratory, 

CWiC denotes the Committee on Machine Methods of Computation, 
DDL denotes the Division of Defense Laboratories, 

100* Comprehe nsive System of Service Routines, developed by the 38EC Group at 
the Digital Computer Laboratory for the input conversion of suitably 
prepared punched paper tapes* When so requested, these routines 
automatically provide a program with suitable programmed arithmetic, 
cycle-counting, and output facilities. 
:DCL Staff: Arden, 12 hours; Best, 31 •75 hours; CombaLic, 6 hours; Demur jian, 
8«5 hours; Denman, 16 hours; Frankovich, 26 hours; Kazel, 16 hours; Helwig, 
80 hours; Kopley, 3 hours; Porter, 21 hours; WWI, 463 minutes 

Testing of the new CS continues. Both the new programs and the old program 
are presently recorded on magnetic tape unit 0, Conversion of- a program by means 
of the new CS can be obtained only if it is specifically requested, 

A memorandum (M-2553) has been written describing the differences between 
the new (30-Jtj) PA and the, old. Programmers who plan to use the new PA should 
read this memorandum before using the program, 

Work has been started on a PA post-mortem for the revised (30-j,j) PA, 

A description of the programs in the CS which involve PA is being written 
and will be included in CS manual 3« 

Helwig 

All of the modifications previously described for the "new" CS are now 
working and are included in the "new" CS, In addition, a modification has been 
successfully tested which will enable CS to convert OCTAL programs. The details 
are described in a forthcoming Bulletin Board Memo, One minor change in the 
unassigned flad table remains to be made, however, before work can be begun on 
the "Direct" CS, 

Frankovich 

The Flad Table routine was tested along with J, Frankovich' s test, for 
unassigned flads. The floating address assignments may be printed in either 
the octal or decimal number system. 
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Corrections and improvements were made on the drum post-mortem subroutines* 
Three tests were made using different entry points to the subroutines. 

Hazel 

The new PA worked for (2^,6) numbers in each of two test routines. These 
test routines will be modified to test other number systems, especially (15,15) 
and (30,0). 

Programs are being written to generate every type of PA alarm in order 
to test the PA-PM, 

Best 

The CS introductory programming course was resumed on December 7 with an 
enrollment of 1? students. Of these students, seven are working on problems 
that have already been accepted for solution on VWI, two expect to subnit pro- 
blems in the very near future, two are from Group 61, and the remaining six 
are associated with the S&EC Group# 

The course consists of 20 one-hour lectures given by members of the S&SC 
Group over a period of two X'jeeks, The material for the first half of the course 
has been described in some detail in the eight completed chapters of GS Manual I 
which is serving as a text for the course. It is expected that the remaining 
ten lectures will soon be written up to complete CS Manual I. 

It is planned to repeat the course during the first two weeks of each 
month provided there is a sufficient number of suitably qualified ap|)licants« 

Denraan, Kopley, Porter 

101 G. Optical P roperties of inpi Metal Fi3jTi? on transparent backings are 

determined and printed out automatically by this program; the input 
data consist of the observed reflection and transmission coefficients, 
the index of the backing, the wavelength, and the sample thickness. 
The program calculates by means of an iterative procedure and prints 
out the index of refraction and the absorption coefficient of the 
film, the rate of variation of these constants with reflection and 
transmission, and the film's conductivity and dielectric constant* 
:for Professor L, Harris, Chemistry Department, Dr, A,L,Loeb 
; bv Dr . A . L ^ Loeb ^ ( DIG ) . 15 hours; J, Richmond (DIG), 20 hours; 
:DCL Staff: Denman, 1 hour; WI, 59 niinutes 

Experimental results published by Krantkramer for thicknesses do^m to 
17a were used as input data, and even the thinnest samples gave satisfactoxy 
runs, in contrast to the saraples reported on by Goes* It was discovered that 
Goos' experimental setup was such as to attempt to eliminate multiple reflection 
in the backing. His results for thinnest films were therefore less than the 
minimum results predicted by the standard equation, so that no set of n and k 
could be found to satisfy Goos' experimental values for the thinnest films. 
This is the reason for the apparent lack of conversion of the program for 
thinnest films. The program now runs satisfactorily for the entire range of 
values tested so far. 
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106 C» MIT Seismic Pro.ject is concerned with the development of methods for 
locating deep reflections from underground strata in seismic 
prospecting. The basic method is one of prediction by means of 
an optimum linear operator* 
sfor Professor PoM^Hurley, Geology and Geophysics; Professor G« Wadsworth, 
Mathematics Department • 

sby E Jl « Robins on ( Res. Assoc,), H. Briscoe, 36 hours; S* Simpson, 20 
hours; W. Walsh, 25 hours 
sDCLs WWI, 46 minutes 

In the last two weeks the gpoup has concentrated on programming rather 
than production. The forward solution of symmetric simultaneous equations has 
been successfiilly coded, using Gauss Elimination and (24,6) arithmetic, and is 
being refined. 

Recent emphasis in the group has been on viewing Linear Operators as 
filtering mechanisms. This stuc^ involves new programs using . Fourier rather 
than correlation techniques, and these programs are under way. 

108 C» An Interpretive Program is being developed that will accept algebraic 

equations, differential equations, etc. expressed on Flexowriter 
punched paper tape in ordinary mathematical notation (within certain 
limits imposed by the Flexowriter) as input and automatically provide 
the desired solution. 

sfor Dr. J,H,Laning,Jr., Instrumentation Laboratory 

sbyDr. J.H.Laning, Jr.(DIC), N. Zierler(DIC) 

:DCL Staffs Hazel, 2 hours; Siegel ; WWI, 12 minutes 

Two problems were solved with very satisfactory results using the 
interpretive routine. 

One of these produced the solutions of a number of independent quadratic 
equations, printing both of the real roots or the conjugate imaginary roots as 
the nature of the equation required. This problem was solved as a test of the 
routine. 

The second problem involved the evaluation of a polynomial function, 

P(n)=p (n) (.2L^zf^2r^^ 
^ ^^ (2L+3)(2L+2)^(2L+l)^ 

with Pq('17 ) = 2 y^ and y^ =0»5 i'or several values of L, The evaluation of this 
function is connected with the work being done for Problem 122, 

109 C« An Airplane Pursuit-Course Program is being developed which will take 

account of airplane dynamics and projectile ballistics and thus 
determine an airplane pursuit course in three dimensions. The 
problem consists essentially of solving 14 simultaneous non-linear 
differential equations by the Runge-Kutta Method which is of fourth- 
order accuracy. 

sfor Mr, J,B,Feldman, Instrumentation Laboratory 

I by M.H.Helljnan (DIC). 40 hours 

jDGL Staff s. Porter, 2,5 hours; WfeTI, 99 minutes 
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A horizontal two-dimensional test run requiring 30 minutes and 5^0 
steps was made to compare the accuracy of the txTO- dimensional horizontal pursuit 
course with that of the three-dimensional pursuit course. The results obtained 
indicate the two-dimensional horizontal pursuit course is a good approximation 
of the three-dimensional pursuit course whose initial conditions are in the 
horizontal plane. 

It is planned to make a number of three-dimensional pursuit course 
runs with the airplane and target initially in a slant plane. From the results 
of these runs it will be possible to evaluate the errors in the calibration 
of existing gunsights, 

113 C e A Stress Analysis of an L-shaped Homogeneous Planar Structure ' is being 
made for the case of a concentrated static load. This struc ture 
is approximated by a framework of bars which will deform in the 
same manner as the prototype. This framework is then analyzed 
using the principles of virtual work and Southwell relaxation 
techniques. Boundary conditions have been specified for the edge 
of the framework so that the deformations of the model will conform 
to the actual deformations of the structure. 

:for Professor J, S.Archer, Department of Civil and Sanitary Engineering 

;by S. Sydney (Res, Assist, CMC), 60 hours 

:DCL Staff: Kopley, 2 hours; mi, 359 minutes 

The program has been tested in sections and is operating satisfactorily. 
The individual units will be combined into a master tape and tested as a whole 
before production runs are made. 

122 B. Coulomb Wave Functions . Regular and irregular solutions of the radial 
Schro'dinger equation with a coulomb potential are being sought. 
These solutions must have the proper asymptotic form to correspond 
to scattered waves from such a potential. The solutions are written 
as a power series in the radial variable and the coefficients are 
determined by recursion formulae, 

:for Profes-sor H, Feshbach, Professor P.M,Morse , Physics Department 

;by A, Temkin (Res. Assist,CMC), 1 hour 

:DCL Staff: Siegel, .5 hour; 'S^I, 9 minutes 

A program has been constructed for finding Pj/1 ), a quantity which 
appears in the recursion relation for the power series coefficients a*J, The relation 
for Pl(t^) is 

?jiri) ^ 1 • (l^Tj^X/f'^T;)^) , , , , . (L^^rf) 2^^ 



^ (2L+1) [jeD l] 



The program tested the value ^J =.5 and went from L=0 to L-6, The latter 
value was 9»6xlO"-L^ and since the values of >] ■.^30»the program shows that this term 

can be left out of the recursion formula for values of L^6, It is my immediate 
prospect to write a program to find the arj . 
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131 o The Training of New Personnel. Tours and Demonstrations are among those 

activities included in this problem. Generally speaking, any approved 
staff problem relating to training and/or demonstrations is considered 
to be in this category. 
sDCL Staffs Kopley . 2 hours; WWI, 37 minutes 

Seven students from Professor Crandall's class at MIT joined Professor 
Cranddll in a tour of the WWI installation on December 1 between 5 and 6 P«M» 
Four other interested outsiders joined the tour which included computer and 
Fl^xowriter demonstrations, 

132 C* Subroutines for the Numerically Controlled Milling Machine are being 

revised and tested. The set of subroutines facilitates programming 
of the computations involved in the preparation of numerical data 
used to control the milling machine. The subroutines involve 
routine numerical and logical operations, 

J for J«0,McDonough, Servomechanisms Laboratory, DIG Project 6873 

•by J,H,R-imyon (E,E,Res. Assist,), 10 hours 

jDGLs WWI, 25 minutes 

A routine for finding points on asymmetrical series I6 airfoil cross- 
sections was successfully tested. Also successfully tested was a routine for 
finding offsets normal to plane curves. This process is used in the first- 
mentioned routine. 

Routines for cut-length selection for constant tolerance and for finding 
the first two derivatives of the series I6 section are still being tested, 

13^ G, Numerical Diagonalization Procedure. This program computes the eigen- 
values and eigenvectors of a symmetric matrix by a method of successive 
rotations. The program is available for use in any problem in which 
this calculation is required, - 

:for Professor J,C .Slater, Physics Department 

I by A, Meckler( DIG). 1 hour 

sDCL. Staffs. Arden, 12 hours; WWI, 90 minutfes 

The secular equation routine is being modified so that the main program 
will be independent of particular data and can therefore be converted once and for 
all. The matrices will be on a separate tape and will be fed in automatically 
for production runs. 

Also, 165 5 :i^ 5 matrices are being diagonalized for Professor J, G, Slater 
and Dr. G.F.Koster, who are doing an electron energy band calculation for body- 
centered structures, 

136, Matrix Equations > ¥arious methods have been, studied for the solution 

of a set of linear algebraic equations, A variation of the Hestenes- 
Stiefel conjugate gradient method has been programmed and tested 
for insertion in the S&EC Library of Subroutines. 
sDCL Staffs Arden, 1^,5 hours; WWI, 81 minutes 

E, Graig of the Electrical Engineering Department is conducting tests 
to detennine the effects of bad matrix donditioning on the accuracy of the 
gradient method that was developed into a subroutine as part of this problem. 
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An attempt is being made to solve a system of order 30 nsing this 
subroutine , 

137 Do investigation of Atmospheric Turbulence as a noise input to airborne 
control systems. A stationary random process is assumed so that 
the methods of generalized harmonic analysis may be used to describe 
the turbulence components in terms of their power spectral densities, 

;for Professor RoG«Seamans, Department of Aeronautical Engineering 

J by Ro A. Summers (Res, Assijfet.)» 6 hours 

sDCLs ma, l4l minutes 

About one-half of the required Fourier transfonm runs have been com- 
pleted and it is anticipated that about one more hour of computer time will 
complete problem 137- 

140« Summer Session System consists of a conversion program, an interpretive 
routine, and mistake diagnostic routines stored in WWI. A special 
mnemonic instruction code has been cfevel oped for use with this system 
thus simulating a computer with characteristics quite different from 
those of W¥I. This Summer Session (SS) computer's developed for 
the use of students participating in the MIT 1953 summer session 
course on "Digital Computers and Their Applications" « The SS computer 
is being used in the E.E*Department courses 6.535 and 6*25 and is 
available to programmers with suitable problems, 
sDCL Staff: Bestj, 33»5 hours; Combelic, 26 hours; Frankovich, 2,5 hours; 
Hoy, 2,5 hours; Siegel, 47 hours; WX, 516 minutes 

Considerable effort has been expended to eliminate errors in the SS 
conversion program and in the instruction rin . both of which must convert numbers 
which appear as a sequence of Flexowriter characters on a punched paper tape into 
a form suitable for storage in an SS memory location. The reasons for the failure 
of these programs to treat numbers containing "xlO^" have been determined, and 
the routines will be corrected in the near future. In addition, a correction has 
been made in the conversion program to permit the use of *'a7i3l "» etc., at the 
end of an SS program to modify the contents of a register which had been filled 
incorrectly in the body of the program. Heretofore, only registers which had 
been assigned a floating address or whose absolute address was known could be 
changed directly by using this device. 

Much WWI time was devoted to running SS programs written by students in 
Course 6o25, The operation of the SS computer was entirely satisfactory. For 
the convenience of these students, the SS computer was temporarily modified to 
display automatically a set of calibrated axes and Vthe names of the programmers 
on the oscilloscope and to phQtjdg;caghi these on the same frame which was to contain 
the curve plotted by each student's program. 

The automatic title-display routine has been completed and incorporated 
into the SS computer. The title of the program is now displayed on all output 
media selected by instructions appearing in the program. 
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141 o S&EC Subroutine Study has been undertaken for the final testing of 
subroutines selected for incorporation into the Library of Sub- 
routines. Although very little effort is going into the specific 
development of subroutines, programs that have been written for 
other S&EC problems and seem to be of general use will be suitably 
modified for the S8iEC Library, 
iDCL Staff: Best, 3 hours;; Denman, 21 hours; Siegel, .25 hour; Vanderburgh . 
1 hour; WWI, Ik minutes 

Since the CS computer will soon be altered to process OCTAL as well as 
DECIML programg,the problem of svdldlng. a dual Subrcutine' Library ha§ come up. In 
order that only one copy of a Library Subroutine need be filed, the following 
measures were adopted: 

All Library Subroutines will start with the characters LSR and 
will end with the words END OF SUBROUTINE ♦ All Library Subroutines 
will be written with DECIMAL addresses. The symbol LSR will indicate 
to the conversion routine that whatever follows^ is DECIMAL, The syiabol 
"END OF SUBROUTINE" will indicate that the original number base is 
to be used again. 

Personal subroutines submitted with the request "I^e in Subroutine 
form" will be typed with LSR at the beginning and END OF SUBROUTINE . at 
the end. They will be filed under the programmer's own tape n-omber* 

Library Subroutines will be filed in a separate file under their LSR 
number, Routines under test for inclusion in the LSR will be filed under problem 
3.^1" in the 'Main. File*. It is proposed that LSR routines be filed in eight major 
categories according to the code shown below: Integer numbers will be used, 
within the categories. The objective here lis to keep the system as simple 
as possible and yet allow seme indexing other than alphabetical* 

GT. ,,,,,♦ l^ewriter (Direct Printer) 

QD.,, ,,,*Delayed I^pei^riter (Magnetic Tape) 

GTD. ,,,,,Both Delayed and Direct !I^ewriter 

OS*,,,,,, Scope (and Camera) 

FU. r , , . , , .Function Evaluation 

MA.,,,,, Matrix (Inversion, Eigenvalues, etc) 

DE ,,,,,. .Differential Equations 

SP,,,,,, •Special (For any routine that does not fit in the other categories) 

List of Available Subroutines for OCTAL or DECIMAL use . 

Tape Numbers 

LSR Tl (24,6)Print 3313 

LSR TDl (30-J,j) Delayed or Direct 2757 

LSR Dl Format .and Delayed 3217 

LSR D2 Delayed (24,6) Print 3333 

LSR D3 Delayed Deciiaal Integer (single) 3214. 

LSR SI (24,6) GD Scope Column Layout 3323 

LSR S2 Single Length Scope Decimal Integer 3322 

LSR Fl (24,6) e^ 3328 

LSR F2 « X 3330 

LSR F3 " In X 3331 

LSR F4 « sin x,cos x 3332 

LSR F5 " sinh x,cosh x 3329 

LSR F6 Arc sin x 3281 
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LSR Ml Largest Eigenvalue 

LSR M2 Simultaneous Equations 

LSR SPL Auxiliaiy ^ Buffer 
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3168 
2624 
2767 



144 G» Self-Consistent Molecular Orbitals are the optimum choices of linear 
combinations of atomic orbitals determined through a process des- 
cribed as a self-consistent field approximation. The niamerical 
procedure involves matrix- vector multiplications, vector additions, 
and matrix diagonalization* 

I for Professor ^'•C .Slater, Physics Department 

».by Dr. A, Meckler (DDL) , 

sDGL Staff: Arden, 4 hours; WWI, 5 minutes 

The use of different starting points indicated no difficulty with 
convergence of the iteration scheme. The program is therefore considered 
suitable for general use. However, it is planned to modify it so that the 
desired data can go on a separate tape and the main program can be converted 
once and for all independent of any particular pre-set parameters, 

147 C, Energy Bands in Crystals are being studied by finding solutions of the 
corresponding second order linear differential equation satisfying 
boundary conditions at the origin. The solutions are found 
approximately by using the Gauss-Jackson formula for forward 
integration. The solutions and their first derivatives are to be 
combined in a sum, the weighting factors being functions of an 
independent parameter. 

sfor Professor J, G. Slater, Physics Department, DIG Ko, 6853 

;by Dr, D,J.Howarth (DIG). 8 hours 

sDGL Staffs Arden, 3 hours; WWI, 535 minutes 

Production is proceeding at a most encouraging rate. The function 
FCEjEq) has been computed for E = -1,8 to +5*0 at intervals of 0*2, and for 
small values of Eq. The roots of this function are required to give a functional 
dependence of Eq on E. For most E, the roots correspond to small Eq, which have 
now beencfetermined. For certain values of E, however, the curve of Eq against 
E has an asymptote. Hence around these E, the calculation will be carried out 
using a smaller range of E and using larger Eq, The results so far obtained 
confirm theoretical predictions of the dependence of Eq upon E, 

1^ C, D1{[1tfl1 Methods of Detecting Signal from Noise are being investigated, 
A sequence of binary numbers will simulate the message wherein 
regions of high density of ones are signal regions and those with 
low density of ones are noise regions. Various methods of detecting 
the change from one region to another, as well as the length and 
midpoint of the signal regions are being studied, 

J for J, f ♦Harrington, Lincoln Laboratory 

' ^J G«P3inneen . (DDL) , 5 hours 

sDCL Staffs Denman, 2 hours; HWI, 23 minutes 

A print-out of the distribution of ones in the binary input demonstrated 
that the pseudo-random-number generator is working satisfactorily. A density 
detector was tried which actually counts the number of ones in a sixteen bit 
interval. This proved. to be the best detector and will serve as a -criterion 
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of judgement « Emphasis is now being placed on designing a detector which 
changes its logic when it detects a target. 

152 D« Diffusion in an Oxide-Coated Cathode is being studied to determine the 

effectfe'of combined thermal and el^cti^Olytic diffusion .that occur in an 
oxide-coated cathode when current is caused to flow through the 
cathode • 

sfor Wo B, Nottingham, Physics Department, DIG No, 63^5 

•^y H.B,Frost (Res. Assist, E.E.Dept,), 12 hours 

jDCL Staff: Porter, 2 hours; WWI, 273 minutes 

During the early part of this period tape 3^09 was debugged completely. 
In this last period three one-hour runs have been completed. A partial run 
of |- hour will soon be completed. To make it easier to carry on partial runs, 
a routine has been incorporated into the program to punch out partial results 
in 55^ form for re-insertion at a later date* 

The 3^9 form of this problem is much faster than the 33II form, mainly 
due to an improved square root routine which makes use of a better first 
approximation. The improvement in speed is by a factor of about 6. 

After each hour run, results differing from the known steady-state 
values by about 3 parts per 1000 were obtained. One run will be extended 
to check the convergence more exactly by comparing the steady-state value 
determined by the program with the known 'value. The truncation error is about 
5 parts in 10,000 depending upon the parameter being run« 

153 G . Gust Response of a Flexible Swept-Wing Airplane is to be "detenained 

for various values of wing loading functions, aircraft configuration 
and dynamic condition parameters^ as input data, giving dynamic 
output data determining the effect of wing flexibility on gust 
response. The solution involves the calculation of forcing functions 
and the evaluation of Duhamel integrals by numerical methods. 
Approximately 120 pairs of linear integro-differential equations 
are to be solved, 
sfor Professor T,H,H,Pian, Aeronautical Engineering Department , DIG No. 6691 
sby K. Foss(DIC),80 hours; D...,Sternlight(Hath.Dept) , 50 hours 
iWU WWI, 58 minutes 

The design specifications of aircraft require a knowledge of the 
acceleration to which an airplane is subjected when flying through a gust(updraft) ; 
consequently^ there has been much active research on this problem. The present 
problem will make a parametric study of the gust response of flexible, swept- 
wing airplanes using a method j^reviotisly developed. 

The two simultaneous integro-differential equations which are to be 
solved are ^ © 



(continued on next page) 
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^A'^(<^)I%o-^)cl^ = ex / f'c^) Y, C^o-^)olcs' 



(1) 



^. fl t^. i'-^ ^;.%, ^j%i u^ i^^^<'!>o-^:> A ^ -*■ A; f^ ^ ^ '*^^'°-^^ ^^ 

^ o (2) 

where the prime mark denotes differentiation with respect to s , the dimension- 
less time-distance parameter. Since solutions are desired only for qj and q^ 
(particularly ql ) the solution will not involve numerical dif f erention j 

only numerical integration. 

The independent parameters of Eqs. (l) and (2) are^ » y^ « » a-., and 
a2s aQ is a constant. The I^ ' functions are known functions of s and {^ j 
where S is an additional independent parameter. The "^ functions -are known 
functions of s^^*^^. and^^ , where ^ is an independent parameter. The F'C^J*) 
function is a known function of Sq and s_, where s-. is an independent parameter*, 

The Duhamel Integrals which occur in Eqs, (l) and (2) have the general 
forms So 

o 
If the interval between and s is divided into m equal stations of width €z ^ 
so that Sq == m-€ ^ the above integral can be evaluated by forming the products 
of z and x at successive stations and integrating them by the trapezoidal rule. 
Then 






O 

From the initial conditions of the problem 



b^r^^ '=' 
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ThuSs the Duhamel integrals can be written in matrix notation for successive 



values of Sq as 



2C<J')x(so-<^^c*-'^^^ 



^'^o 



■iJai \t 



O 

o 



■«h-i 



X, 



m-A. 



X 



rh^3 



► ♦ IT' 



o 
o 

o 

o 



t-'^<. 



■m 



L 



J 



Bqs^ (1) and (2), therefore, can be written in matrix notation as 



(5) 



(6) 



(7) 



where each of the matrices [k], [B], [cj .and f DJ is of the same general form 
as Eq» (5) J specifically, they are triangular matrices, each column of which is 
foimed by shifting the preceding column down one row (only the first column 
need by computed to form the matrix) • 



For the problem at hand, there will be 20 intervals of numerical integration^ 
each of length s /lO^ which yields matrices of order 20 x 20, Because of the 
triangular matrices involved, Eqs, (6) and (7) can be developed into relatively 
simple recurrence formulas for each value of q" or q9 , in terms of all previous 
values • There will b e one or two thousand such solutions for different combinations 
of the independent parameters, /, ^ , s >/^ ,/^ , a^ ^and a^* 

For each different combination of/, fi and s , the functions I V, \^ . 
and FK^) must be recomputed^ The function F'(^) is simply 






(8) 
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The I^ ' function is given by 









The other I functions are similar, each containing the same logarithmic term. 
The expansion of this logarithmic term about unity into a Taylor's series results 
in 






^L^[v^ [^±^^\^£4J^ f" -J — IJtiS (10) 



The termination of this series at r = 4 has been found to give sufficient 
accuracy for all the values of f and d" which will be used. 

The ^ functions are similar but more complicated functions, and they 
must also be expanded into a series for numerical evaluation. 



During this perio,d programs. w^i^Sr^'Wi'itten ^. a tested for evaluating the 
Ip n^ 'H' « and f functions. The V and f furidtions were , successfully evaluated 
after certain parameters had been suitably scale factored. The program for 
evaluating the I functions Jias not given satisfactory results. Further study 
is being carrieci out.' '*' 



where 
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155 B* Synoptic Climatology « A multiple regression fornrala is used to predict 
temperatures fran. pressure distributions described by Tschebycheff 
polynomials* The matrix of scalar products which is used in the 
calculation of the coefficients of the multiple~regression system is 
being calculated on WSI.. 

:for Professor T.F»Malone, Meteorology Department 

*by R. Miller (DIC) 

sDCL Staff; Arden, 13 hours; Porter, .5 hour; WWI, 26 minutes 

A (30,0) decimal integer print routine has been successfully corapleted 
for use in the program. The data for the first production run is being prepared* 

157 C» Rectangular Matrix Multiplication is being programmed for inclusion in 
the S&EC Library of Subroutines* 
sfor Dr. A, Meckler, Solid State Physics 
sby Dr. A. Meckler(DDl,), 4 hours 
sDCLs WI, 23 minutes 

The multiplication of two rectangular matrices has been successfully . 
programmed and the routine is now included in the S&EC Library of Subroutines. 
The main program using this routine must plant the various orders (which should 
not exceed 32 ) of the matrices and the drum groups on which the matrices are 
stored. The multiplication routine then carries out the calculation in, (2^,6) 
leaving each matrix factor on its drt^m group and the product on another. 

1*3 Operating Statistics 

1.31 Computer Tme 

The following indicates the distribution of Wl time allocated to the 
S&EC Group. 



Programs 
Conversion 
Magnetic Drum Test 
Magnetic Tape Test 
Scope Calibration 
Demons trations( #131) 
Total Time Used ;■ 
Total TiRie Assigned' - 
Isable Tirrie:^ Percentage 
Number of Programs 



53 hours, 
13 hours, 



Id 

10 

19 
26 
52 



minutes 
minutes 
minutes 
minutes 
minutes 
minutes 



68 hours, 34 
74 hours, 12 
92.2$^ 
238 



minutes 
minutes 



1.32 Program Time Distribution 

The following table attempts to show how the WSI time expended on S&EC 
programs was distributed with respect to machine runs that gave meaningful results 
(productive computer time) and runs that gave unsatisfactory results (♦♦lost" 
computer time}. Productive computer time is subdivided to indicate the time 
involved in actual computations as contrasted with the time expended getting 
information out of WWI. Computer time "lost" is subdivided to show the portion 
of time lost due to errors in the programmer' s f omiulation of his problem 
(logical errors); due to errors in the programmer's use of the WWI code, CS 
Conventions, etc. (technical errors); due to tape preparation errors; due to 
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errors by the S&EC computer operators in running the program; due to malfunctioning 
of terminal equipment; and finally due to miscellaneous causes* 

These times are determined as percentages of the time listed.' above in 
section 1#31 for programs. The figures below have been averaged over the last 
three biweekly periods. The times used in computing these figures are extracted 
from the biweekly report forms subnitted by the various programmers who have 
used S&EC allocated WWl time, 

1. Productive Computer Time 
Computation 55*^ 
Output 9.9/^ 

2. Computer Time Lost Due to Programmers* Errors 
Technical 19*5$ 

Logical 4.;^ 

3. Computer Time Lost Due to Other Difficulties 
Tape Preparation 2«7^ 
Operator's Errors 2.0$S 
Terminal Equipment Malfunction 2,2^ 
Miscellaneous 3*6^ 

1*33 Tape Preparation (M, Mackey) 

An attempt is being made to obtain some idea of the time expended in the 
preparation of tapes* During the past biweekly period a check was made on the 
tapes processed* 

Due to the variations in procedures involved we have distinguished kmbng 
original complete tapes and the following three types: typed modifications - 
changes of 11 or more registers which must be typed, converted, then attached 
to the main program or changes which must be made in the body of a Flexowriter 
tape; manual modifications - changes punched directly?- in 556 form and attached 
to a converted tape; combined tapes - which require duplication of two or more 
complete tapes. 

The following information was compiled: 

Complete Typed Manual Combined 

Tapes Mods Mods Tapes 

No. of Tapes l4l 55 48 79 

No. of Registers 38,213 1277 2^7 

Time Consumed 115 hrs.35 min. 27 hrs.53 min. 8 hrs.ll min, 10 hrs. 10 min, 

Thus, it may be seen that the average length of an original complete tape is 
271 registers requiring 49 minutes to prepare. A typed modification averages 
22 registers in length and requires 29 minutes to prepare while Manual Modifi- 
cations average 6 registers and require 10 minutes for preparation. 
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2. COMPUTER ENGIHEERING 

2.1 WWI System Operation 

2.11 Core Memory (L.D.Holmes, A.J.Roberts) 

The drum parity-check system has been installed and is now functioning. 
Drum parity alarms should occur only on the bi or rd orders. Core Memory alarms 
should not occur on these orders. 

2.12 Auxiliary-Prum System (K.E,Mc Vicar) 

The parity system has been installed for the auxiliary drum, and parity 
checks are now made on all transfers from the drum. 

An attempt is being made to revise the paperwork for the drum so that it 
fits into the Whirlwind system. This includes the compilation of parts lists, 
making of wiring schedules, and use of WWI modification notices for all drum- 
system changes. 

2.13 Magnetic Tape (E.P.Farnsworth) 

Minnesota Mining and Manufacturing has finally received a shipment of mylar 
from Dupont and expects that 80 percent of the tape we have on order will be 
shipped splice-free. They have just developed a new splicing tape having less 
than one-tenth the stretch of the standard kl splicing tape; splices made with 
a sample of the new tape appear to be usable but may cause trouble after a 
few months storage or use, Three-M has also advised that tri -acetate base tape 
is not adequate for our application. Audio- tape is shipping us two sample reels 
of unspliced high-output mylar-base tape. 

The adapters necessary for converting Unit 2 to triangular-^hub reels were 
completed by the shop and are ready for installation. Panels for operating a 
second printout on Unit 2 are being installed. Replacement of a pair of electro- 
lytic capacitors in the servo-chassis which caused instability with Unit 
previously and Unit 2 presently has apparently cleared the trouble. These cap- 
acitors are being checked for intermittent leakage* 

A trouble-report sheet has been prepared, and the marginal clutch test 
program has been modified to correct for the change in computer speed. 

2.14 Typewriter and Paper Tape (L.H.Korcott) 

Recent carriage-return troubles on the delayed- output printers were apparently 
caused by extra ** start tape" pulses originating in the magnetic-tape printout- 
control register. Farnsworth has modified his circuits to correct this condition. 

2.2 Terminal Equipment 

2.21 Ferranti PETR (F.E.Irlsh) 

Sketches for the mechanical and electronic assembly of the amplifier cir- 
cuits of the Ferranti PETR and a sketch of the control circuits have been given 
to the drafting room. 
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According to the present schedule, one Ferranti PE1R will be installed 
in a temporary fashion during the next two weeks. The constiniction of the 
final amplifier units for the three Ferranti PElR's and the control circuits 
for one unit will be completed by February 1, 1953* The final installation 
will follow that date, 

2.3 General 

2.31 WWI Service File (D.A. Morrison) 

Plans are under way to relocate the WWI Service File in WrTL Test Control 
to allow the Flexowriter console to be moved out of the passageway. The 
contents of the steel cabinets are to be consolidated into a smaller cabinet, 
and the Wffl Sejrvice File and the smaller cabinet will occupy the space now 
used fcy the steel cabinets. The upper left-hand drawer of the desk in Test 
Control now contains the test tapes. 



3» LIBRARY ACCESSIONS LIST 

The following material has been received in the Library, ¥2-325 • 
Library Files 
No. Source Title 



2592 Business Week, Reprint Computers Go Commercial 

2593 Ji®S Methods of Conjugate Gradients for Solving 

Linear Systems 

2598 Rensselaer Computer Lab. A Magnetic Tape Recorder for Analog Computing 



^- ADMINISTRATION AND PERSONNEL 
Terminated Staff ( J. C. Proctor) 
Richard Farmer 

Terminated Non-Staff (R. A. Osborne) 

Assunta Aprile 
Nancy Jones 
Athena Korologos 
Esther Sidman 



